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Abstract. We investigate the newly discovered supersolid phase by solving in random-phase approximation 
the anisotropic Heisenberg model of the hard-core boson 4 He lattice at zero temperature. We include nearest 
and next-nearest neighbor interactions and calculate exactly all pair correlation functions in a cumulant 
decoupling scheme. We demonstrate the importance of vacancies and interstitials in the formation of the 
supersolid phase. The supersolid phase is characterized by strong quantum fluctuations which are taken 
into account rigorously. Furthermore we confirm that the superfiuid to supersolid transition is triggered 
by a collapsing roton minimum however is stable against spontaneously induced superflow, i.e. vortex 
creation. 



PACS. 05.30.Jp Boson systems - 67.80. -s Quantum solids 
4 He - 75.10.Jm Quantized spin models 



67.80.bd Superfluidity in solid He, supersolid 



1 Introduction 

The counterintuitive idea of a superflow in a solid, later 
coined supersolidity was firstly conjectured in 1969 by 
AndreevpQ and in 1970 seized by Leggett and Chester[21 
[3]. From a theoretical point of view, supersolidity is a 
state of matter characterized by simultaneous off-diagonal 
(ODLRO) and diagonal long range order (DLRO). It was 
speculated that such a phase exists because vacancies and 
other defectons are non-localized and will Bose condense 
at sufficiently low temperature. Still most physicists re- 
mained critical of the notion as several experiments failed 
to produce any evidence of this state. Finally in 2004 Kim 
and Chan 4, 5 measured a tiny superflow in solid he- 
lium at temperatures below T=0.2 K, expressed by non- 
classical rotational inertia in a torsional oscillation exper- 
iment, and thus proved the existence of the supersolid 
state. This landmark experiment rekindled vast interest 
in the supersolid state and subsequently many new the- 
ories and numerical quantum Monte-Carlo calculations 
supporting the existence of supersolidity were proposed. 
However the true nature of the supersolid phase still re- 
mains obscure. Numerous follow-up experiments managed 
to shed light on the matter but the relevance of 3 He im- 
purities and especially the nature of the unconventional 
normal solid to supersolid transition resembling the 2D 
Kosterlitz-Thouless transition is still being debated. 
Recent experiments [S] raised new questions as it was found 
that the supersolid phase exhibits a hysteresis, where the 
superfiuid signal depends on the chronology of variation 
of temperature and in the amplitude of the rotational os- 
cillation. An other recent experiment detected a change 



in the elastic properties of solid helium. The change of the 
elastic moduli bears a remarkable resemblance with the 
supersolid signal. 

However, despite sophisticated numerical methods and ad- 
vanced theories such as vortex liquids [5] and superglass 
states [5] we believe that there remains a gap in the range 
of theories of the supersolid phase. In this paper we in- 
tend to fill this gap and present a theory of supersolidity 
in a quantum lattice gas (QLG) model beyond classical 
mean-field. We follow the approach of K.-S. Liu and M.E. 
Fisher [10] and map the QLG model to the anisotropic 
Heisenberg model. The method of Green's functions proved 
to be very successful in the description of ferromagnetic 
and anti-ferromagnetic states and we use this method to 
investigate the supersolid phase which corresponds to a 
canted anti-ferromagnetic phase. 

The emergent third order Green's functions in the randor 
phase approximation (RPA) are broken down using the 
cumulant decoupling to yield a closed set of equations. 
Quantum fluctuations at zero temperature result in va- 
cancies and interstitials present even at zero temperature 
and in the supersolid phase the net vacancy density is 
therefore non zero. The supersolid phase is characterized 
by Bose condensation of the vacancies as well as the inter- 
stitials and thus both will contribute to superfluidity. In- 
terestingly, the major contribution comes from vacancies. 
Also, our model confirms propositions that the superfiuid 
to supersolid transition is triggered by a collapsing roton 
minimum j 1 1 [1 1 2 j . Nonetheless our solution shows, contrary 
to earlier results that this transition is stable against spon- 
taneously induced superflow. 

The paper is organized as follows: In Section [2] and [3] we 
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introduce the generic Hamiltonian of a bosonic many body 
system and discretize it to a quantum lattice gas model. 
This model is equivalent to the anisotropic Heisenberg 
model in an external field and we will identify the cor- 
responding phases. In the following two sections we derive 
basic thermodynamic properties relevant at zero tempera- 
ture and discuss the excitation spectrum of the spin waves 
in the superfluid and the supersolid phases. The relevance 
of quantum fluctuations is discussed in Section[7l where we 
provide justification of results briefly reported elsewhere 
[13 . Finally in the last two sections we discuss key prop- 
erties of the supersolid phase and their occurrence within 
the example of three different sets of coupling constants. 



2 Generic Hamiltonian 

Apart from possible 3 He impurities the supersolid Helium- 
4 state is a bosonic system and the generic Hamiltonian 
for such systems in the language of second quantization is 
given by: 



H 



d 3 xipUx)(-^-V 2 + fiWx) 
2m 



d xd x'iJj* (x)ip* (x')V(x — x')ip(x)'ip(x') 



(1) 



where tp^(x), the particle creation operator and ip(x) , the 
corresponding destruction operator obey the usual bosonic 
commutator relations. Hamiltonians in three dimensions 
such as in Eq. ([1]) are not solvable even for elementary po- 
tentials V(x) such as the Dirac delta distribution. There- 
fore we are induced to introduce further approximations. 
An approximation which proved particularly successful for 
the description of liquid Helium is know as the Quantum 
Lattice Gas model and was first introduced by Matsubara 
and Matsuda[i"4]. 

In the quantum lattice gas model one works with a space 
lattice of discrete lattice points rather than the contin- 
uum. This approximation shows to be very useful for the 
supersolid state as the spatial discretization of the this 
model serves as a natural frame for the crystal lattice of 
(super)-solid helium. This procedure significantly simpli- 
fies the problem of breaking translational invariance sym- 
metry for states that exhibit diagonal long range order. 
In this way this model gives the easiest possible access to 
analyze states that exhibit both diagonal and off diago- 
nal long range order simultaneously. Also in this model 
no specific knowledge of the density distribution of the 
atoms is needed. 

According to Matsubara and TsunetofIS] the generic Hamil- 
tonian Eq. ([I]) in the discrete lattice model reads: 



Here Uij are non-zero for nearest neighbor and next near- 
est neighbor hopping and otherwise zero. The values of 
u nn and Unnn are such that the kinetic energy is isotropic 
up to the 4th order. In the case of a BCC lattice (two in- 
terpenetrating SC lattices) the matrix elements are given 
by: 



2 1 



3 4roa 2 
1 1 



3 4ma 2 



(3) 



As atoms do not penetrate each other there can exist only 
one atom at a time on a lattice site and consequently 
and a are the creation and annihilation operators of a hard 
core boson commuting on different lattice sites: 



(i + 3) 



(4) 



but obey the anti-commutator relations on identical sites: 

[4> a t\+ = i a i> a i\+ = ° 



di, a. 



1 



(5) 



Equation ([2]) is the Hubbard model in 3 dimensions for 
hard core bosons. Due to the unusual statistics of hard 
core bosons, there does not exist a Wick's theorem for 
their operators and the common formalism of pertubative 
field theory is not applicable. Hence in the following chap- 
ter we transform the model to an equivalent spin model 
namely the anisotropic Heisenberg model. 



3 Anisotropic Heisenberg Model 

It is well known that the operators of hard-core bosons 
obey the same SU(2) algebra as spin S — 1/2 particles 
do. Therefore it is feasible to replace the creation and an- 
nihilation operators by spin operators. 



2 j - s j 
2 j = 



y - 2 s j- 



(6) 



It is easily apprehensible that the usual lie algebra for 
spin 1/2 particles preserves the mixed commutation/anti- 
commutation relations for hard-core bosons. This substi- 
tution transforms the hard-core bosonic Hubbard model 
into a spin model: 



i 



2SfS? - 2SfSf) 



si 

2 



St 



(7) 
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Fig. 1. The BCC lattice consists of two interpenetrating SC 
sub-lattices. In the perfect solid phase one sub-lattice (i.e. sub- 
lattice A) serves as a lattice of on-site centers and is occupied 
while sub-lattice B represents the empty interstitial and is va- 
cant. 

if we adjust the notation to conform with the usual stan- 
dards of spin models, it becomes evident that the resulting 
Hamiltonian is the anisotropic Heisenberg model : 

i 

with: 



ft* = A*-E J «+E J Jl (9) 

3 3 



Table 1. Possible magnetic phases and the corresponding 
phases of the Hubbard model. All Phases are defined by their 
long range order. The columns, from left to right, are the spin 
configurations, magnetic phases, ODLRO, DLRO and corre- 
sponding 4 He phases, respectively. 



TT 


FE 


No 


No 


Normal Liquid 


// 


CFE 


Yes 


No 


Superfluid 


// 


CAF 


Yes 


Yes 


Supersolid 


U 


AF 


No 


Yes 


Normal Solid 



ferromagnetic and anti-ferromagnetic phases which we ab- 
breviate by FE, CFE, CAF and AF. The order parame- 
ters, mi for off-diagonal long range order and mi for diag- 
onal long range order in the magnetic system are defined 
by: 

TOl = (51) + (S B ) 

m 2 = (S A ) - (S' B ) (10) 

In the following we will use these order parameters to iden- 
tify the phases within the classical mean-field approxima- 
tion as was derived by Matsuda and Tsuneto [15] and 
extended to finite temperature by Fisher and Liu [10] as 
well as in the novel random-phase approximation. 



4 Green's Function 



The Hubbard model as well as the anisotropic Heisenberg 
model are defined on a discrete lattice and one may ask 
to what extend a specific choice of lattice geometry will 
affect the physical properties of the system. While the 
quantitative results certainly depend on the lattice geom- 
etry we can safely assume that qualitative properties, such 
as phase transitions and critical constants will not change 
for different lattices as long as no frustration effects are 
evoked. Therefore we may safely choose, in order to avoid 
unnecessary complications, a simple lattice geometry and 
an obvious choice are two interpenetrating simple cubic 
sub-lattices which together form a BCC lattice, see Fig- 
ure [TJ Defining two sub-lattices gives us the possibility 
to establish the DLRO of solids in a natural way: sub- 
lattice A represents the centers of the 4 He ions, hence it 
coincides with the ion lattice. Sub-lattice B defines the 
interstitial, the space in-between those atomic centers. In 
the liquid phases of course the occupation number on both 
sub-lattices is equal as there is no spatial density variation. 
Table [T] charts the various magnetic phases and identi- 
fies the corresponding phases of the 4 He system. Accord- 
ing to the spin configurations we call the four magnetic 
phases ferromagnetic, canted ferromagnetic, canted anti- 



The anisotropic Heisenberg model in an external field is, 
not least due to absence of 0(3) symmetry difficult to 
solve. However, in the context of supersolidity, investi- 
gating the CAF phase, it has been solved in a classical 
mean-field approximation 10, 15j. The anisotropic Heisen- 
berg Hamiltonian in the classical mean-field approxima- 
tion is obtained by substituting the spin-1/2 operators 
with their respective expectation values: 

H MF = ~h z ((S A ) + (S' B )) 
-2jf(S A )(S B ) - 4((S A )(S A ) + (S B )(S B )) 
-2J?(S%)(S%) - Jj((S A )(S A ) + (S%)S%)), (11) 

where jj 1 = qij\ eAjeB , j{ = Qajf^Aj&A > J i = <?i JJeAjeB 
and Jj = <liJj e Aj£Ai 9i = 6 and qi = 8 are the number of 
nearest and next nearest neighbors. Although the classical 
mean-field theory is quite insightful and gives an accurate 
description of the variously ordered phases its fails to take 
quantum fluctuations and quasi-particle excitations into 
account. Hence, in order to overcome these shortcomings 
we derive a fully quantum mechanical approximation and 
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CFE 



CFE-CAF 
transition 



CAF 



AF 



Magnetic Field h z 

Fig. 2. The possible mean-field phases of the anisotropic 
Heisenberg model on a bipartite lattice with external field h z . 



solve the anisotropic Heisenberg model in the random- 
phase approximation which is based on the Green's func- 
tion technique. At zero temperature the retarded and ad- 
vanced Tyablikov [TTUTE] commutator Green's function 
defined in real time are: 



&-(*) 
(*) 



ijAd 



-i6(t)(no\[Snt),S%}\no) 
ie(-t)(n \[S?(t),S%]\n ) 



(12) 



Here \no) is the normalized Heisenberg ground state, /i 
and v are elements of {x, y, z} and i and j denote the lat- 
tices sites. The most successful technique of solving many 
body Green's function involves the method of equation of 
motion which is given by: 



id ^LM = s(t)([s?,s?])-ie(t)([[sf,H) 



idtG 



XI) 



(t) = 6(t)([s?,sy]) + ie(-t)([[s? 



(13) 



The commutator [Sf , H] can be eliminated by using the 
Heisenberg equation of motion giving rise to higher, third 
order Green's functions on the RHS. In order to obtain a 
closed set of equations we apply the cumulant decoupling 
procedure and as a consequence the third order Green's 
functions split into product terms of single operator ex- 
pectation values and second order Green's functions. The 
cumulant decoupling |20j is based on the assumption that 
the last term of the following equality is negligible: 

(Abc) = 

(A)(BC) + (B)(Ad) 

+ (C)(Ab)- 2(A) (B)(6) 

+ ((A-(A))(B-(B))(C-(C))) (14) 

This decoupling scheme leads to a closed set of six equa- 
tions that determine the six Green's functions, correspond- 
ing to three spin components on two sub-lattices. 
The cumulant decoupling scheme introduces mean-fields 
of the spins operators which have to be calculated in a self- 
consistent manner. The Green's functions gives us the pos- 
sibility to calculate a set of two self-consistency equations. 
However, in the supersolid and supcrfluid case the off- 
diagonal long range order; i.e. non-zero transversal fields 
(S A ) and (S A ) increase the degrees of freedom in number 
by two and therefore analytical properties of the commu- 
tator Green's functions pose two additional constraints on 



the spin-fields given by: 



h* + 2(S%)(4 - Jj) + 2(S B )j\ = 23} 



h* + 2(S Z B )(4 - Jj) + 2{S* A ) J{ = 2jj 



(s%) 
(si) 

(S X A ) 



(S%) 



(S z b) (15) 



These equations are quite important and also hold in the 
classical mean-field approximation, where the state of the 
system is defined by these equations together with 
V(S X A ) 2 + (S A ) 2 = \ and v/(S!> 2 + <S|) 2 = \. Equation 
(fl~5|) also determines possible second order phase transi- 
tions as was shown my Matsuda and TsunetofiT)]. The 
normal fluid to superfluid (FE-CFE) second order transi- 
tion is located at: 



n FE-CFE — J l 



J 2 jj' jj 



(16) 



For the normal solid to supersolid (CAF-AF) the classical 
mean-field approximation yields: 



L CAF-AF 



fi- 



4+4 -m 2 -(4 



(17) 



and the critical magnetic field h z (corresponds to the chem- 
ical potential /i in the language of the QLG) for the su- 
pcrfluid to supersolid (CFE-CAF) transition is: 



l CFE-CAF 



J I jj Ji + J J 

y/ (-4 + 4 - JJ) 2 - (J?r- 



(18) 



For a particular choice of coupling constants all four phases 
will exists when, see Figure [2] 



l FE—CFE 



>h 



CFE-CAF 



>h z r 



CAF-AF 



(19) 



The three equations for the critical fields, Eq. (fT6"|) . Eq. 
(|17p and Eq. (fl"8|) , are derived for in classical mean-field 
approximation. However as Equation ()15|) also holds in 
the random-phase approximation these values give a good 
indication where the actual transitions take place. Never- 
theless, due to depletion of the spin-fields caused by quan- 
tum fluctuations the actual values are slightly lower, see 
Section El 



5 Thermodynamic Properties 

In the first section we have seen that the grand-canonical 
QLG Hamiltonian, where the number of particles are vari- 
able corresponds to the canonical anisotropic Heisenberg 
Hamiltonian with fixed number of spins. Therefore there 
exists following relation between any thermodynamic po- 
tential defined in the QLG model and the anisotropic 
Heisenberg model: 



Oqlg - = H 



eisenberg 



(20) 



Here refers to an arbitrary thermodynamic potential. 
The self-consistency equations, derived in the previous 
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sections, determine the spin fields of the anisotropic Heisen- 
berg model in the various phases, but in regions of h z 
where more than one solution may exist, we have to com- 
pare internal energies to select the true ground state. In- 
tuitively, we might want to compute the internal energy 
by calculating the expectation value of the Hamiltonian 
as: 



H=h*j2( s i)+J2 J U s i s j) 



(SIS')) 



(21) 



where the correlation functions can be derived from the 
corresponding Green's functions, Eq. (fT2"|) . Unfortunately 
this approach will yield incorrect and inconsistent results 
as the Green's function derived with the cumulant decou- 
pling is not an exact solution of the anisotropic Heisenberg 
Hamiltonian but rather the solutions of an unknown, effec- 
tive model, which is an approximation of the Heisenberg 
model. As we do not know the exact form of this effective 
model we have to resort to fundamental thermodynamic 
relations to integrate the energy. At T=0 following equal- 
ity holds: 

au _ (s z A ) + (s* B ) 

dh z 2 



(22) 




I )| 



[100] 



[110] 



[111] 



I ] 



This equation allows us to determine if the critical fields 
given by Eq. ([15]). Eq. (T7]) and Eq. (H]) really refer to 
second order phase transitions. As the z-component of the 
spin is decreased in the canted ferromagnetic phase due 
to the onset of the transversal field this phase is energeti- 
cally favorable over the ferromagnetic phase. Similarly the 
canted anti-ferromagnetic phase has lower energy than the 
anti-ferromagnetic phase as the total magnetization in z- 
direction in the canted anti-ferromagnetic phase is some- 
what higher due to the increasing influence of the ferro- 
magnetic term h z J7 Sf . Hence, the phases as depicted in 
Figure [2] are real. 

Unfortunately Eq. (|2"2"|) only allows one to determine the 
energetically favorable state when the possible transition 
point is known, such as in second order phase transitions. 
But we can not use the relation to determine possible first 
order transitions and unfortunately this shortcoming is 
only resolvable by extending the formalism to finite T. 
Thermodynamic relations connect various macroscopic quan- 
tities and we will use them to obtain observable properties. 
Although the external magnetic field h z in the spin model 
is an observable the corresponding quantity in the QLG 
model, namely the chemical potential is not. Therefore we 
are interested to attain a formula for the pressure associ- 
ated with a certain chemical potential. The relationships 
is most easily derived from the following Maxwell relation: 



Fig. 3. Excitation spectrum of the anisotropic Heisenberg 
model with = 1.498^, jj = 0.562if, jf = -3.899K and 
jj = — 1.782_K" for the supersolid phase (solid line and dashed 
line) at h z = 0.65 and the superfluid phase (long dashed line) 
at h z = 7.46. Here [000], [100], [110], [110] refer to the various 
points of the first Brillouin zone. 

6 Excitation spectrum 

The excitation spectrum in the ferromagnetic phase as 
well as in the anti-ferromagnetic phase feature the well 
known, gaped magnon excitation with quadratic k depen- 
dence in the k — > limit. In the canted ferromagnetic 
(superfluid) phase the spin-wave excitation spectrum is 
comprised, due to spontaneously broken U(l) symmetry, 
of the gapless linear Goldstone phonons. 
Additionally the canted anti-ferromagnetic (supersolid) phase 
exhibits a second branch which accounts for the broken 
translational symmetry. This second branch goes quadratic 
with k in the long wavelength limit and has a gap: 



A = [4 2(4 - 4 
2 + (S Z A ) 2 



(S x a) 



Jj)(s x A ) 2 (s x B ) 2 

(S x a) 3 



(S x b ) 2 + (S Z B ) 2 



(S X B ) 



(24) 



fdP\ fdN\ #, ottice 

sites \ 
\dH) T ,y = \dv) T ^ = V [ 



(23) 



where e := (S A ) + (S B ). These are the basic thermody- 
namic properties that we will use in the further discussion 
and in principle all other properties can be derived from 
the internal energy U. 



In Figure [3] the quasi-particle excitation spectrum is plot- 
ted for the superfluid (CFE) and the supersolid (CAF) 
phases. In the supersolid (CAF) phase the excitation en- 
ergy reaches zero at the edge (point [100]) of the first Bril- 
louin zone. The corresponding spin-waves refer to oscilla- 
tions with a 7r-phasc shift between different sub-lattices. 
Hence, on a single sub-lattice the spin-wave looks like a 
zero wave-number mode. It was recently suggested |llj 
that the superfluid to supersolid transition is triggered 
by a collapsing roton minimum. This assumption is sup- 
ported by the present model; here the transition to the 
supersolid phase takes place when the excitation spectrum 
goes soft at [100]. The dispersion relation in the superfluid 
(CFE) phase is given by: 



o;(fc) = 2{(J 1 T ( 7l (fc)-l) 
[<S z ) 2 (Ji T (7i(fc)-l) 



•/ 2 T (72(fc)- 
- J 2 T ( 72 (fc) 



l))x 
1))- 
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(S X ) 2 (JJ + Jj - Jhi(k) - Jh(k))] } 1/2 , (25) 

where 71 (k) and 72 (k) denote the standard lattice gener- 
ating functions of a BCC lattice. From this equation we 
can see that the energy possibly goes to zero at [100] (cor- 
responds to 71 (k) — — 1 and 72 (fc) = 1) when the following 
condition is fulfilled: 



j -I + j 2 ' + J{ - 4 1 < 



(26) 



Hence we obtained a further condition (supplementary to 
Eq. (TT8|) ) for the existence of the superfluid to supersolid 
transition. 

Equation (f25f allows for the existence of a second region 
of the reciprocal space where the dispersion relation might 
go soft. For 7i(fe) = and 72(fc) = —1 which corresponds 
to [111] we obtain following condition: 



-2J 2 T 



J J + 



> 



7" 

Jo 



It was also conjectured that the superfluid phase in the 
vicinity of the superfluid to supersolid transition is un- 
stable against spontaneously induced superflow and su- 
perflow associated with vortices. Therefore we investigate 
how the present model reacts to induced superflow. A net 
superflow is either given by a moving condensate which 
results in a gradient of the phase of the wave-function, or 
equivalently by a moving environment while the conden- 
sate stays at rest. The latter is obtained by an additional 
term in the Hamiltonian: 



d 3 x<if> 1 (x)(iv n ■ V)ip(x) 



(28) 



which corresponds to following term in the anisotropic 
Heisenberg model: 



ill 



^J*(S?SV-SfS?) 



(29) 



Here the nearest and next-nearest neighbor cross coupling 
constant are anti-symmetric Jy = —Jfi and are zero for 
directions perpendicular to the motion of the environ- 
ment v n . The term yields an additional matrix M± in he 
random-phase approximation: 



M 1 = (J?{k) + J}(k)) 



(S z ) ' 
(S z ) 
-(S x ) 



(30) 



The matrix is reduced to dimension 3x3 because we are 
only interested in the superfluid phase only where no dif- 
ference between the two sub-lattices is made. The cross 
coupling terms J x are given by: 



J* = ^ J ab exp(ifcaAs) 

aAB 

J* = ^ JAA ex P( ika AA), 



(31) 



where ciaa and aAB are the lattice parameters correspond- 
ing to A and B lattice sites. The dispersion relations, given 
by the eigenvalues of the total matrix M + M\ , are altered 
in the following way: 



uo k -> (J 1 x (/c) + J*(k))+LU k 



(32) 



in the k — > limit this accounts for a tilt of the dis- 
persion curve toward the motion of the environment; the 
quasi-particle energy in the direction of the motion v n is 
lowered while the energy for particles traveling in oppo- 
site directions is lifted. From the definition of the lattice 
generating functions Ji(k) and 72(A) and Eq. (|31|) we see 
that J* (A) = J 2 x (/c) = for k where 71 (k) = -1 and 
72(A) = 1. Hence the roton dip that triggers the superfluid 
to supersolid transition is not affected by the superflow. 
The situation is likewise for the roton minimum at [111] 
where 71 (k) = and 72 (k) = — 1. Also, here the cross 
couplings J* and J 2 X become zero and the stability is not 



(27) affected by induced superflow. 



7 Quantum Fluctuations at Zero 
Temperature 

As presented in Ref. [13] quantum fluctuations are impor- 
tant to study as they can lead to the understanding of 
the physical origin of the different phases. At zero tem- 
perature there are no thermal fluctuations present in the 
system and all fluctuations will stem from quantum me- 
chanical effects. The mean-field approximation as derived 
in the beginning of the paper is a classical approxima- 
tion and as such it does not display quantum fluctuations. 
This is expressed by a constant spin magnitude of 1/2 
over all phases at zero temperature. In the anisotropic 
Heisenberg model quantum fluctuations are a result of 
non- vanishing pair correlations (5f 5J) — (Sf ) (Sj) of near- 
est and next nearest neighbors. Consequently, as random- 
phase approximation takes those correlations accurately 
into account we expect quantum fluctuations which are 
expressed in a depletion of the total spin magnitude as 
can be seen in Figure 2J In the ferromagnetic phase the 
total spin is 1/2 and thus there are no quantum fluctua- 
tions present. This is expected as the ferromagnetic phase 
is governed in the h z — > 00 limit by an effective single 
operator, and hence the classical Hamiltonian: 



H 



(33) 



The spin depletion is strongly pronounced in the CFE 
and CAF phases where transversal components account 
for additional fluctuations. We also see that in the CAF 
phase the spin magnitude is different on the two sub- 
lattices. This indicates that next nearest neighbor inter- 
actions are dominant and balancing nearest neighbor in- 
tegrations are slightly suppressed. Therefore we assume 
that in the canted anti-ferromagnetic phase the two sub- 
lattices do slightly decouple. This interpretation is sup- 
ported by the spin- wave excitation spectrum in the canted 
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2 4 6 8 

h 1 

Fig. 4. The total spin magnitude at zero temperature is shown 
for all four phases, for the cases J± = 1.498A", jj = 0.562A', 
jf = -3.899A and jJj = -1.782A, main plot and jj = 0.5K, 
Jj = 0.5A, j| = -2K and j| = -0.5A, inset, respectively. 
The mean-field solution gives always the horizontal dotted line. 
In RPA the total spin is depleted due to quantum fluctuations. 
For the supersolid (CAF) phase the RPA gives a difference 
between the sublattices A and B (solid and dashed lines main 
plot). For details, see text. 

anti-ferromagnetic phase. We have seen that there exists 
a zero frequency mode, where the spins on different sub- 
lattices are 7r-phase shifted. Hence, spin-waves which cou- 
ple both sub-lattices carry no energy. 

8 Discussion 
8.1 Case 1 

In this section we will discuss the properties of the super- 
solid phase using the example of two sets of coupling con- 
stants. As we are interested in describing real systems we 
may ask what sets of parameters are physical and which 
set exhibits the best fit to 4 He. Physically, the transver- 
sal constants J T ought to be positive as they correspond 
to the kinetic energy. In quantum lattice gas models usu- 
ally J T s are chosen so that the kinetic energy is isotropic 
up to 4th order giving the best possible approximation 
to the continuum limit. However in the supersolid phase 
the Hamiltonian may be regarded as an effective model 
and therefore we refrain from this restriction. The inter- 
actions between the helium atoms are controlled by van- 
der-Waals forces and their repulsive nature at very short 
distances determines negative nearest neighbor interac- 
tion j| , evoking anti-ferromagnetic ordering in the spin 
language. The corresponding Lennard- Jones potential is 
short ranged and therefore it is sufficient to only consider 
nearest and next nearest neighbor interactions. Liu and 
Fisher chose coupling constants in order to fit the model 
to the actual phase diagram of Helium-4. As the supersolid 
phase had not been discovered experimentally at this time 
they investigated the possibility of a supersolid phase ex- 
isting. 




H 1 U U U.J 1 l.J 



Fig. 5. Figure a) shows the relation between pressure in the 
QLG model and the external magnetic field in the anisotropic 
Heisenberg model. The pressure is normalized with respect to 
the critical pressure P c marking the supersolid to superfluid 
transition. The dashed line refers to set 1: = 1.498A, jj = 
0.562A, jj 1 = -3.899A and jj = -1.782A and the solid line 
to set 2: jj = 0.5K, Jj = 0.5K, jj 1 = -2.0K and jJj = 
— 0.5A Figure b) depicts the internal energy of the model as a 
function of the pressure for set 1 (dashed line) and set 2 (solid 
line) . 

It is widely accepted that the lambda transitions falls into 
the universality class of the XY-model, which is a limiting 
case of the anisotropic Heisenberg model. In the same way 
we believe that the anisotropic Heisenberg model is capa- 
ble of covering the essential properties of the supersolid 
phase. Nevertheless the present model will fail to appro- 
priately map 4 He over the complete range of temperature 
and pressure as various properties such as variability of 
the lattice constant and lattice vibration modes (phonons) 
are not taken into account in this model. Therefore we ab- 
stain from fitting the solutions of random-phase approxi- 
mation to the phase diagram of real 4 He but merely choose 
two sample sets to study key properties of the supersolid 
phase. 

The first set of parameters is given by: 

Jj = 0.5K 
J 2 T = 0.5K 
j\ =-2K 

jj = -.5K (34) 

In the classical mean-field this set of parameters exhibits 
all four phases where the corresponding critical magnetic 
fields are given by Eq. (JIBJ, Eq. IJT5J) and Eq. {T7]) and 
yield respectively: h z FE _ CFE — 3.5, h z JFE _ CAF — 2.0207 
and h z CAF _ AF — 0.86608. In the classical mean-field as 
well as in the random-phase approximation the transitions 
are determined by Equation (|I5p . which shows that h z 
roughly scales with S A and S%. Therefore we expect that 
the transitions in the random-phase approximation due 
to depletion of the spin magnitude are slightly lower. The 
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actual values are: h z FF _ ( j FF = 3.5, h z cFE—CAF = 

1.96 

and h,Q AF _ AF = 0.857. 

The second set we have chosen was extensively studied by 
Liu and Fisher and their parameters are given by: 



Ji = 


IA98K 


Jj = 


0.562X 


4 = - 


•3.899X 


4 = - 


1.782 A" 



k03 



(35) 

Again the transition points are slightly lower (except for 
the FE-CFE transition) for the random-phase approxi- 



mation and yield: h z F _ CFE 
1.0071(1.0577) and h z CAF _ AF 



7.741(7.741), h 



CFE-CAF 



0.3963(0.41716), where 



).85 0.9 0.95 



the numbers in parentheses are the values derived by the 
classical mean-field approximation. Figure [5] a) depicts 
the relation between the external field of the anisotropic 
Hcisenberg model and the pressure in the QLG as given 
by Equation (|23p . The pressure on the y-axis is renormal- 
ized such that 1 corresponds to the critical pressure P c of 
the superfluid to supersolid transition, given by roughly 20 
atm in Helium 4. High magnetic field corresponds to low 
pressure and vice versa. Negative magnetic fields corre- 
spond to high pressures that do not have physical validity 
in the quantum lattice gas. Therefore the maximal pres- 
sure corresponds to zero magnetic fields. 
In Figure[5]b) we plotted the internal energy of the anisotropic 
Heisenberg model which corresponds to the conjugated 
potential C/[/i] = U — fiN per volume of the QLG model 
which is minimized at zero temperature. In Figure [HI we 
plotted the superfluid order parameter (ip) as a function 
of the pressure in the superfluid and supersolid phase for 
both sets of parameters. The order parameter displays its 
maximum value in the vicinity of the transition to the su- 
perfluid phase and evidently approaches zero at the NS-SS 
transition. The superfluid order parameter on the on-site 
sub-lattice associated with vacancies is higher than the one 
on the interstitial sub-lattice. While this effect is strongly 
pronounced in Set 1, where the order parameter of the 
vacancies is around 37 times the order parameter of the 
interstitials near the SS-NS transition, in Set 2 the Bose 
condensation of the vacancies is only marginally higher 
(1.3 times) than of the interstitials. Yet we observe that 
in this model Bose condensation appears in the vacancies 
as well as in the interstitials though the major contribu- 
tion comes from the vacancies. In Figure we have the 
density of vacancies, the interstitials and the difference of 
both, the net vacancy density plotted as a function of the 
pressure in the supersolid and normal solid phases. We see 
that in the normal solid phase the number of vacancies and 
interstitials stays finite. This is due to quantum fluctua- 
tions and consequently the number of vacancies have in 
addition to thermal activation a second contribution re- 
sulting from quantum mechanical effects. Nevertheless the 
number of vacancies and interstitials appear in equal num- 
bers and the net contribution in the normal solid is zero. 
This is different in the supersolid phase where a surplus 
of vacancies accounts for a positive net vacancy density. 
As we decrease the pressure in the supersolid phase both 



1 1.05 

p/p. 



Fig. 6. 



Magnitude of the superfluid order parameter 
the supersolid and superfluid phase. Long dashed line refers 
to the on-site sub-lattice A while the dashed line refers to the 
interstitial sub-lattice B. The dotted line is the average of both. 
Figure a) refers to set 1 and Figure b) to set 2 (see text). 




Fig. 7. Density of vacancies (dashed line), density of inter- 
stitials (long dashed line) and the difference of both, the net 
vacancy density (solid line) in the normal solid (for P/P c > 1) 
and the supersolid phases (below 1). Figure a) refers to set 1 
and Figure b) to set 2 (see text). 



the vacancy density and the interstitial density increase. 
However, the vacancy density increases faster, leaving a 
net vacancy density which reaches its maximum at the 
phases transition between the supersolid and the super- 
fluid. Interestingly the net vacancy density varies nearly 
linearly with the pressure as the solid line in Figure [7] 
shows. 



8.2 Case 2 

In the section on the excitation spectrum we have seen 
that the spin-wave energy at [100] of the first Brillouin 
zone goes soft exactly when the superfluity to supersolid 
phase transition occurs. Additionally, for coupling con- 
stants that fulfill condition Eq. (12711 there exists a second 
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and we will leave the detailed discussion of this phase to 
future work. 




[000] [100] [110] [111] [000J 

Fig. 8. Excitation spectrum (solid line) in the superfluid phase 
for Jj = 0.5K, Jj = 0.5K, jf = -2K j| = -1.5K just before 
the phase become unstable due to a collapsing minimum at 
[111]. Induced superflow alters the spectrum (dashed line) but 
does not affect the minimum. 



minimum at [111] which can collapse. Following set of con- 
stants fulfill this condition: 



Ji = 


0.5K 


Jj = 


0.5K 


J 11 - 


-2K 


4 = - 


-1.5K 


n (IHJl 


there 



superfluid) transition in the system: 



FE—CFE 



4.5 



(36) 



(37) 



beneath this line classical mean-field approximation pre- 
dicts a CFE (superfluid) phase that extends down to h z = 

as due to the relatively large negative jj the solid phase 
does not acquire a sufficiently low free energy to be the 
true ground state. The random-phase approximation how- 
ever draws a slightly different picture. Analogous to the 
classical mean-field solution the random-phase approxima- 
tion also yields a phase transition near h z — 4.5. But un- 
like the classical mean-field solution, the superfluid phase 
here does not survive all the way down to h z = 0. Due to 
the particular choice of parameters the superfluid phase 
becomes unstable at around h z = 2; i.e the quasi-particle 
spectrum turns imaginary at 7i(fc) = and 72 (k) = — 1 
([111]), as Figure [8] shows. The dashed line in this figure 
shows the excitation spectrum under an induced super- 
flow v s . The roton minimum is not affected by this super- 
flow and hence the superfluid phase is not destabilized by 
an spontaneously induced superflow. Interestingly beyond 
this line no other stable phase exists in the random-phase 
approximation; there is no set of spin fields (SJ), (Sg), 
(S%) and (Sg) that solves the self-consistency equations of 
the random-phase approximation. Therefore we conclude 
that there must exist a 'novel' phase that is not covered 
by the random-phase approximation on a bipartite lattice 



9 Conclusion 

In this paper we analysed the supersolid phase in the three 
dimensional quantum lattice gas model. Through transfor- 
mation to the anisotropic Heisenberg model in a external 
field we were able to employ the well-established tech- 
nique of real-time Green's functions for spin systems. The 
series of infinite order Green's functions as it appears in 
the equation of motion was truncated by applying cumu- 
lant decoupling and the resulting random-phase approxi- 
mation accounts for linear spin- waves. We are the first to 
apply this method to the canted anti-ferromagnetic phase 
entailing a set of 6 algebraic equations. The innate self- 
consistency equations inhere a 3 dimensional numerical 
integral over the k-space. By introducing a two dimen- 
sional density of states the integral was reduced to two 
dimensions where the lattice generating functions serve as 
new integration variables. In the said integral the DOS is 
the only quantity that depends on the structure of the lat- 
tice. Hence, once the DOS is computed for a certain lattice 
geometry the further calculation remain unaltered. There- 
fore our method is widely applicable and easily adjustable 
to various magnetic systems where canted phases are in 
the center of interest. This also holds for 2 dimensional 
lattices where linear spin waves are expected to yield a 
reasonable approximation. 

The random-phase approximation takes quantum fluctu- 
ations into account and consequently in this solution the 
solid phase exhibits vacancies and interstitials at zero tem- 
perature. Yet in the normal solid phase the vacancies and 
the interstitials occur in equal number, thus yielding a zero 
net vacancy density. In the supersolid phase this balance 
shifts in favor of the vacancies giving rise to a finite pos- 
itive net vacancy density at zero temperature. Our data 
also shows that vacancies as well as interstitials Bose con- 
dense and hence both contribute to superfluidity. Never- 
theless the Bose condensation is stronger expressed in the 
vacancies thus giving the major contribution to superso- 
lidity. 

Furthermore the present approach confirms suggestions 
that the superfluid to supersolid transition is triggered by 
a collapsing roton minimum. However our results show 
that this roton dip is not affected by Galilean transforma- 
tion and hence the superfluid phase is stable against spon- 
taneously induced superflow. Additionally we find that for 
a narrow regime of parameters a second roton minimum 
collapses. Below this point a stable phase does not exist 
in the bipartite random-phase approximation and a so- 
lution is thus beyond the model. The prospect of future 
work looks promising. The formalism is easily extendable 
to finite temperatures as shown in Ref. ^22], where we in- 
vestigated the properties of the supersolid phase at finite 
T. In particular the temperature dependence of the net va- 
cancy density and the behavior of the specific heat across 
the supersolid to normal solid transitions is of particular 
interest. 
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